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Cosmological perturbations with wavelengths smaller than Hubble 
radius can be handled in the context of Newtonian theory with very 
high accuracy. The application of this Newtonian approximation, how- 
ever, is restricted to nonrelativistic matter and cannot be used for rela- 
tivistic matter. Recently, by modifying the continuity equation, Lima, 
et. al., extended the domain of applicability of Newtonian cosmology 
to radiation dominated phase. We adopted this continuity equation 
to re-examine linear cosmological perturbation theory for a two fluid 
universe with uniform pressure. We study the evolution equations for 
density contrasts and their validity in different epochs and on scales 

qq , larger than Hubble radius and compare the results with the full rel- 

■ ativistic approach. The comparison shows the high accuracy of this 

^\ ' approximation. 
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I. INTRODUCTION 



The real universe contains inhomogeneous structures like galaxies, clusters etc and, in any theory of the formation 
of these structures, it is essential to understand the evolution of small inhomogeneities in the early universe. In 
' principle, it is straightforward to work out the general relativistic theory of linear perturbations j^,^). By linearizing 
Einstein's equations, we can obtain a second-order differential equation of the form 

t> ; 

£(9ap) $9at3 = ST a fi, (1) 

00 , 

where C is a linear differential operator depending on the background metric, g a f3, and the set (Sg a 0,ST a ^) denotes 
CJ the perturbations in the metric and stress tensor about an expanding FRW background. 

In practice, however, there are many complications and coneptual difficulties which make this analysis highly 
, nontrivial. One issue is the so-called "gauge problem" (extensively discussed in the literature; see [3-13]) which arises 
60- due to nonuniqueness of splitting the metric and matter variables into a zeroth order background and small, first-order 
perturbations. Since a relabeling of coordinates x a — > x a can make a small 5T a p large or even generate a component 
which was originally absent, one must care to factor out effects due to coordinate transformations, when analyzing 
relativistic perturbations. There are two different ways of handling these difficulties in general relativity. One approach 
is to analyze a perturbation in a particular gauge, say, synchronous gauge fl~4|| . In this case, one specifically identifies 
the points of a fictitious background spacetime with those of real spacetime, and treat ST°q to be the perturbed 
mass density etc. In this method, however, we cannot fix the gauge completely and the residual gauge ambiguities 
can create some problems. The second method is to construct the perturbed physical variables in a gauge-invariant 
manner. The gauge-invariant approach is conceptually more attractive since there is no need for specific identification 
of the points between the two spacetimes, though it is more complicated and the physical meaning of variables do 
not, in general, posses any simple interpretation and becomes obvious only for specific observers. 

It is convenient to divide the cosmological perturbations into two subclasses: (i) Perturbations with wavelengths 
larger than Hubble scale (A > du), for which we have to use some form of a general relativistic theory of perturbations 
and (ii) small-scale perturbations (A < du) for which the evolution of mass density can be studied using Newtonian 
theory. In this context, all physical quantities can be defined unambiguously to the order of accuracy needed. In 
general, the application of Newtonian equations is further restricted to nonrelativistic matter and cannot be used for 
relativistic component even for scales much smaller than Hubble radius (A <C dj{)- 
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Recently, Lima [ |1 5[ , et al. , re-examined the basic equations describing a Newtonian universe with uniform pressure 
and found a way for obtaining the same evolution equation for density contrast as could be obtained by the full 
relativistic approach. They achieved this goal by modifying the continuity equation in an expanding background. 
Using this result, they argued that one can extend the domain of validity of Newtonian cosmology |16| in order to 
analyze some problems of formation of structures even in the radiation dominated phase. 

In this paper we extend the result of reference [flEf to a multi-component universe with different equations of state. 
We shall then consider a two fluid universe in the context of "pseudo Newtonian" cosmology. Comparison with the fully 
relativistic two-component universe reveals the high accuracy of density contrast equations in the pseudo Newtonian 
cosmology. 

The organization of this paper is as follows: In section II we will use the modified continuity equation in order to 
find out the evolution equation for density contrast of a multi-component universe with arbitrary equations of state. 
This is actually the generalization of paper by Lima |l5| , et. al. In section III we will examine the result of section II 
for the two fluid system with radiation and dark matter. In this section we will try to give approximate solutions to 
coupled equations of density contrast of radiation and dark matter. We shall compare the result in this context with 
those ones as obtained in fully relativistic approach, for example given in ref [0. We give the conclusions in section 
IV. 



II. THE LINEAR PERTURBATION THEORY FOR A MULTI-COMPONENT UNIVERSE 

The set of equations describing a several component universe, in the proper coordinates (i, r) with r = a(t)x, can 
be written as, (i) "Continuity" equation, (i\)Euler equation, and (iii) Poisson equation, 

^ + V r • (pu) + pV r • u = 0, (2) 

^ + (u-V r )u = -V,.7/>-(p + j 3 )- 1 V r p, (3) 

V r 2 iP = 4irG(p + 3p), (4) 

where the "continuity" equation is the one that has been used by Lima ]l5| , et. al., for an expanding universe, to 
study the evolution of cosmological perturbation using linear regime. Here, the velocity of light, c, assumed to be 
unity and the quantities 

N N N N 

p = ^2pi, p = ^2pi, u = X] ui and ^ = X]^ ( 5 ) 

i— 1 i— 1 i— 1 i—l 

denote the total density, total pressure, total field velocity, and the total generalized gravitational potential of the 
cosmic fluid. 

In an expanding Newtonian universe the evolution of small fluctuations can be studied in the usual perturbation 
theory, i.e, by perturbing the background density, pressure, field velocity, and gravitational potential. That is, we 
take 

p(r,t)=p bg (t)[l + S{r,t)}, (6) 
p(r,t) = Pbg (t) + 6p(r,t), (7) 
il>(T,t)=il> bg (T,t)+<p(r,t), (8) 
u(r,t) = u 6a +v(r,t), (9) 

where the second term is a small correction. By inserting the above expressions into Eqs (1), linearizing the resulting 
equations to first order in perturbations, and transforming to comoving coordinates by x = (r/a), we get, 

i -i- iA x ) 
#*> + i±r_ V*-vW=0 

a 

*W + *vW = -l V .y- f V,J<*> (10) 
a a (1 + j/Wja 



VbV = 47rGa 2 ^(1 + 3^)pb Si ^ 
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where we have assumed that the expansion is adiabatic and thus vi X ^ = (5p/5p) = v^ x \ is the sound velocity. Here 
(X) denotes the species under consideration and pi = vipi. The sum on the right hand side is over all components. 
Note that v to t and v s 2 ot are not constant any more and thus the relation v s 2 ot = v to t does not hold for the full system. 

By eliminating the peculiar velocity between the above equations, and Fourier transforming the perturbation such 
that V x 2 <5 (x) = -k 2 5 {x) etc., we get 

k X) + > + (^^) 2 & W = 4*G(1 + J*>) + ^) Pbgi 8 H . (11) 

This is a peculiar equation in the sense that there is no counterpart for this equation in the fully relativistic 
treatment, unless we impose the synchronous gauge and comoving gauge simultaneously, (one cannot impose the two 
gauge conditions simultaneously even in large scale limit in the presence of pressure, since we get 5 = system when we 
ignore the entropic and anisotropic pressures [JlSj] .) In the case of single component medium and (k 2 v s 2 /a 2 ) <C 1, 
however, equation ( |il| ) leads to the same equation derived in full relativistic approach by imposing the two gauge 
conditions simultaneously (see for example, eq. (15.10.57) in ref. E3, eq. (10.118) in ref. jl9| and eq. (6.136) in 
ref. M). 



III. TWO FLUID UNIVERSE 



Now, we shall assume that our universe consists of radiation [R] and dark matter [DM1 with vr = 1/3 = v 2 R and 



vdm = "m 



= 0. Then the evolution of 5^(t) and 5j? M (t) can be determined from eq. (JTl|), to give 



a 



5n, + 2-5r 



a 



k 2 v\ 32 



ttG ( Sp DM \ 

~ pr v + ^ ) 



Sr = 0. (12) 



W + 2-5dm - iirGpDM (l + P^-) 
a V OPDM J 



5 dm = 0. (13) 



Though these equations cannot be solved exactly in closed form all the important properties of the solutions can 
be obtained by suitable approximations. We shall now discuss these properties by introducing the new variable 
x = (a/a eq ) where a eq is the expansion factor at the time t = t eq . Transforming from t to x as independent variable 
we can re-express all the quantities in terms of x, 

?™=-L; ?2- = -^; e** = ± i [x + l). (14) 

Peq 2x^ Peq 2x Peq 2x^ 

Since (jptot/ Peq) = {vdm +Pr)/ Peq — (l/6x ) we get the following result for the full two component system 

Vtot Vdm+Vr 1 / 1K s 

v = Vtot = = : = ^77- — 7> ( 15 ) 

Ptot PDM + PR 3(1 + X) 



and 



1 



H 2 (x)=H 2 —(l + x). (16) 



eq 2x 4 



On the other hand, one can express the combination (k / H a ) as 

k 2 2x 2 



(17) 



H 2 a 2 (1 + x) 



where w is defined by the ratio 2ituj = [dn (a eq ) / ^(a eq )}. In terms of this variable, x, the time derivatives can be given 
by 

H{x)x4- = Ht> 
dx 



d 








dt 


da 


d 2 






= H 2 a4- 


dt 2 


da 



^ = H ' 2 ^ I ) " 2^ 2 (1 + v)a Ta = H 2 D 2 - -H 2 (l + u)D, (18) 
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where D — x(d/dx) and we have used the new continuity equation 
equations for Sr and S dm become 



With these modifications, the coupled 



D z 



and 



1 



2(1 + a;) 



D 



2 lo 2 x 2 



3(1 + a;) (l + x) 



Sr = 



2x 



(l + x) 



>DM , 



(19) 



D + 



-D 



2 (l + x) 2(1 



TDM 



(l+x) 



(20) 



A particular mode, labeled by the parameter lo, enters the Hubble radius in the radiation dominated phase if 
lo > 1 and in matter dominated phase if a; < 1. The lo and x en t are related to each other through the relation 
Lo 2 x 2 nt = 2-n 2 (l + x en t). We shall consider the case lo 2 ^> 1 which has more features and is more relevant to this 
paper. Then three cases need to be discussed, (i) When xlo <C 1, in which case the wave length of the mode is bigger 
than the Hubble radius; (ii) when lo^ 1 « i « 1, the mode has entered the Hubble radius in the radiation dominated 
phase; and finally the case, (iii) when lox 3> 1, the mode is inside the Hubble radius in a matter dominated phase. 
Let us now consider each case. 

Case (i) xuj <C 1 

In this case the equation (ft9|) and (Efl) give 



The growing solution is 



D 2 -4)5 R = 0; 



D 2 S DM = 3<5_r. 



(21) 



8dm = x 2 



(22) 



Thus both radiation and dark matter grow as x 2 oc a 2 . Obviously, since the universe is radiation dominated 6 dm is 
driven by Sr and its back reaction on Sr is negligible. This result matches with the fully relativistic one! 

Case (ii) lox ^S> 1 and x -C 1 



The concerned equations in this approximation are 



D z + -wV 



Sr = 0; D 2 5dm 



3S R 



The solution for Sr equation would be of the form 



C\ Jo 



+ C 2 N Q 



(23) 



(24) 



where Jo an d A^o are the Bessel functions of the first and second kind, respectively. If we use the fact that lox 1 
the approximate solution to <5r becomes 



IT 1 " UX , 



(25) 



which is oscillating rapidly since lox 3> 1. Substituting for Sr in the second equation and using the fact that x -C 1, 
we find that 



'DM 



6C 



sm 




A In: 



B. 



(26) 



Therefore the perturbation in dark matter grows essentially logarithmically during the period lo~ x « i « 1, which 
also known from previous work. Though the wavelength of the dark matter perturbation is bigger than effective Jeans 
length, the growth of dark matter perturbation is prevented by the rapid expansion of the universe. The similarities of 
these results with their fully relativistic ones come from using the modified continuity equation, which-in turn-permits 
one to enlarge the domain of applicability of Newtonian cosmology even to radiation dominated phase. 
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Case (iii) lux ^> 1 

This range corresponds to the matter dominated phase with mode inside the Hubble radius. The coupled equations 
now become 



D 



Id- 3 - 

2 2 



Sdm — 0; 



D 



1 - 2 o 

-D + -U3 X 



25dm 



Solving for Sdm we get 



Sdm = Ax + Bx 2/3 = Ax. 
Plugging this into the Sr equation we find that 



dx 



dSj- 



IdS 



2 dx 



3A 



(27) 



(28) 



(29) 



Choosing, Sr = (3A/lu 2 ) + x 3 > i f(x), we get the following differential equation for f{x) 

By applying the WKB approximation, we may solve the / equation and get the final result for Sr to be 



(30) 



3A 



B- 



exp ±i\ 



,1/2 



(31) 



The oscillations presented by the second term, in Sr, continue to dominate over the driving by the Sdm term. Hence, 
the radiation density does not grow in the matter dominated universe, i.e., when Xr < du- These solutions are very 
similar to the solutions to the fully relativistic equations. The reasons are clear during this phase: First of all, we are 
dealing with the modes which are well within the Hubble radius and secondly the modes are in the matter dominated 
phase. 



IV. CONCLUSIONS 



We have shown that the result of the linear cosmological perturbation of a two fluid universe can be obtained with 
high accuracy from Newtonian cosmology even in the presence of pressure and for scales larger than Hubble radius. 
This is obtained by using the modified Newtonian equations. In fact, by using the modified continuity equation one 
can get rid of a misleading pressure gradient term which is obtained in the semi-classical formulation and more over 
obtain the time evolution of the density contrast for any value of parameter v. Comparison of © and © with 
fully relativistic equations, say, in the comoving gauge (see e.g. rcf. p7|) given by: 



D 



2(1 



-D 



2 tu 2 x 2 



3(1 



(1 



5r- 



2.T 



(1 + x) 



Sdm — F(Sr,Sdm), 



(32) 



and 



D 



1 x 



D - 



3 x 



2{l + x) 21 + x 



Sdm - 1 — x 5 R = G(S R ), 



(33) 



with 



F(S Rl S 



dm 



D 



4 6x 2 + 13a; + 8 
9 (z + 4/3) 2 (l + x) 



- --, — X , DSdm ■ 



3 (x + 4/3) 



(34) 



and 



G(S R ) = 



z + 4/3 



A _ 1 6x 2 + 13x 



3(l + x)(rc + 4/3) 



(35) 
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will allow us to estimate the accuracy. Left hand sides of Eqs. ([32]) and ( p3| ) are the same as their pseudo Newtonian 
counterparts Eqs. (|ll| ) and (^0|), while F and G are the corrections coming from general relativity. Numerical 
comparison of F and G with their left hand sides in Eqs. ([l9]) and (|2(]) reveals that these corrections are very small 
and almost negligible indeed. This can be better seen in Figs. 1 and 2. In each of these figures the vertical axis is the 
ratio of pseudo Newtonian density contrast to its exact relativistic value at the same epoch. Both the figures show 
that this ratio for the radiation component is almost unity and the value of dark matter density contrast in pseudo 
Newtonian case differs at most by a factor 2. As we discussed earlier, the solutions of Eqs. (|r]) and (|o|) for three 
different cases are in a very good agreement with the solutions of Eqs. ( |32"|) and (|33|). This agreement for ui ^> 1 is 
more profound (see Fig. 1 and Fig. 2). In other words, the ratio R = (5 Newton/^ exact) f° r radiation tends to unity 
for all epochs as we increase the value of to, though remains off for dark matter by a lesser factor. It seems that the 
Friedmann's equations have strong correspondence with Newtonian theory, even more than one might have naively 
expected. 

AN wishes to thank Robert Brandenberger for many useful discussions. 
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Figure captions 



Fig. 1. Comparison of perturbation amplitudes in dark matter (DM) in pseudo Newtonian cosmology (PNC) and 
relativistic cosmology (RC). The amplitude of DM in NC is more than RC with an error of less than factor 2 for 
uo = 100. After z = 390 up to present time the ratio becomes constant. The radiation amplitude in RC is closer to 
that in NC. The behavior of two perturbations are almost the same except for some small phase difference. The ratio 
R = Newton Inexact) for radiation is almost one except for some tiny fluctuations. 

Fig. 2 This is same as Fig. 1 with w = 200. Note that the agreement between the modes are more pronounced. The 
ratio for dark matter in this case is less than 2 and it becomes a constant after z = 39 . In radiation dominated phase 
when modes are inside the Hubble radius, modes are almost in phase and again amplitude of radiation is more in RC 
compare to PNC at any instant. 
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